Revealing the atomic and electronic mechanism of human manganese superoxide dismutase product inhibition.

Human manganese superoxide dismutase (MnSOD) is a crucial oxidoreductase that maintains the vitality of mitochondria by converting O2∙- to O2 and H2O2 with proton-coupled electron transfers (PCETs). Since changes in mitochondrial H2O2 concentrations are capable of stimulating apoptotic signaling pathways, human MnSOD has evolutionarily gained the ability to be highly inhibited by its own product, H2O2. A separate set of PCETs is thought to regulate product inhibition, though mechanisms of PCETs are typically unknown due to difficulties in detecting the protonation states of specific residues that coincide with the electronic state of the redox center. To shed light on the underlying mechanism, we combined neutron diffraction and X-ray absorption spectroscopy of the product-bound, trivalent, and divalent states to reveal the all-atom structures and electronic configuration of the metal. The data identifies the product-inhibited complex for the first time and a PCET mechanism of inhibition is constructed.

Revealing the atomic and electronic mechanism of human manganese superoxide dismutase product inhibition.

Human manganese superoxide dismutase (MnSOD) is a crucial oxidoreductase that maintains the vitality of mitochondria by converting O 2 ●- to O 2 and H 2 O 2 with proton-coupled electron transfers (PCETs). Since changes in mitochondrial H 2 O 2 concentrations are capable of stimulating apoptotic signaling pathways, human MnSOD has evolutionarily gained the ability to be highly inhibited by its own product, H 2 O 2 . A separate set of PCETs is thought to regulate product inhibition, though mechanisms of PCETs are typically unknown due to difficulties in detecting the protonation states of specific residues that coincide with the electronic state of the redox center. To shed light on the underlying mechanism, we combined neutron diffraction and X-ray absorption spectroscopy of the product-bound, trivalent, and divalent states to reveal the all-atom structures and electronic configuration of the metal. The data identifies the product-inhibited complex for the first time and a PCET mechanism of inhibition is constructed.

Insights into the mechanism of SARS-CoV-2 main protease autocatalytic maturation from model precursors.

A critical step for SARS-CoV-2 assembly and maturation involves the autoactivation of the main protease (MProWT) from precursor polyproteins. Upon expression, a model precursor of MProWT mediates its own release at its termini rapidly to yield a mature dimer. A construct with an E290A mutation within MPro exhibits time dependent autoprocessing of the accumulated precursor at the N-terminal nsp4/nsp5 site followed by the C-terminal nsp5/nsp6 cleavage. In contrast, a precursor containing E290A and R298A mutations (MProM) displays cleavage only at the nsp4/nsp5 site to yield an intermediate monomeric product, which is cleaved at the nsp5/nsp6 site only by MProWT. MProM and the catalytic domain (MPro1-199) fused to the truncated nsp4 region also show time-dependent conversion in vitro to produce MProM and MPro1-199, respectively. The reactions follow first-order kinetics indicating that the nsp4/nsp5 cleavage occurs via an intramolecular mechanism. These results support a mechanism involving an N-terminal intramolecular cleavage leading to an increase in the dimer population and followed by an intermolecular cleavage at the C-terminus. Thus, targeting the predominantly monomeric MPro precursor for inhibition may lead to the identification of potent drugs for treatment.

Perfect Crystals: microgravity capillary counterdiffusion crystallization of human manganese superoxide dismutase for neutron crystallography.

The NASA mission Perfect Crystals used the microgravity environment on the International Space Station (ISS) to grow crystals of human manganese superoxide dismutase (MnSOD)-an oxidoreductase critical for mitochondrial vitality and human health. The mission's overarching aim is to perform neutron protein crystallography (NPC) on MnSOD to directly visualize proton positions and derive a chemical understanding of the concerted proton electron transfers performed by the enzyme. Large crystals that are perfect enough to diffract neutrons to sufficient resolution are essential for NPC. This combination, large and perfect, is hard to achieve on Earth due to gravity-induced convective mixing. Capillary counterdiffusion methods were developed that provided a gradient of conditions for crystal growth along with a built-in time delay that prevented premature crystallization before stowage on the ISS. Here, we report a highly successful and versatile crystallization system to grow a plethora of crystals for high-resolution NPC.

Contribution of the catalytic dyad of SARS-CoV-2 main protease to binding covalent and noncovalent inhibitors.

The effect of mutations of the catalytic dyad residues of SARS-CoV-2 main protease (MProWT) on the thermodynamics of binding of covalent inhibitors comprising nitrile [nirmatrelvir (NMV), NBH2], aldehyde (GC373), and ketone (BBH1) warheads to MPro is examined together with room temperature X-ray crystallography. When lacking the nucleophilic C145, NMV binding is ∼400-fold weaker corresponding to 3.5 kcal/mol and 13.3 °C decrease in free energy (ΔG) and thermal stability (Tm), respectively, relative to MProWT. The H41A mutation results in a 20-fold increase in the dissociation constant (Kd), and 1.7 kcal/mol and 1.4 °C decreases in ΔG and Tm, respectively. Increasing the pH from 7.2 to 8.2 enhances NMV binding to MProH41A, whereas no significant change is observed in binding to MProWT. Structures of the four inhibitor complexes with MPro1-304/C145A show that the active site geometries of the complexes are nearly identical to that of MProWT with the nucleophilic sulfur of C145 positioned to react with the nitrile or the carbonyl carbon. These results support a two-step mechanism for the formation of the covalent complex involving an initial non-covalent binding followed by a nucleophilic attack by the thiolate anion of C145 on the warhead carbon. Noncovalent inhibitor ensitrelvir (ESV) exhibits a binding affinity to MProWT that is similar to NMV but differs in its thermodynamic signature from NMV. The binding of ESV to MProC145A also results in a significant, but smaller, increase in Kd and decrease in ΔG and Tm, relative to NMV.

AI-Accelerated Design of Targeted Covalent Inhibitors for SARS-CoV-2.

Direct-acting antivirals for the treatment of the COVID-19 pandemic caused by the SARS-CoV-2 virus are needed to complement vaccination efforts. Given the ongoing emergence of new variants, automated experimentation, and active learning based fast workflows for antiviral lead discovery remain critical to our ability to address the pandemic's evolution in a timely manner. While several such pipelines have been introduced to discover candidates with noncovalent interactions with the main protease (Mpro), here we developed a closed-loop artificial intelligence pipeline to design electrophilic warhead-based covalent candidates. This work introduces a deep learning-assisted automated computational workflow to introduce linkers and an electrophilic "warhead" to design covalent candidates and incorporates cutting-edge experimental techniques for validation. Using this process, promising candidates in the library were screened, and several potential hits were identified and tested experimentally using native mass spectrometry and fluorescence resonance energy transfer (FRET)-based screening assays. We identified four chloroacetamide-based covalent inhibitors of Mpro with micromolar affinities (KI of 5.27 µM) using our pipeline. Experimentally resolved binding modes for each compound were determined using room-temperature X-ray crystallography, which is consistent with the predicted poses. The induced conformational changes based on molecular dynamics simulations further suggest that the dynamics may be an important factor to further improve selectivity, thereby effectively lowering KI and reducing toxicity. These results demonstrate the utility of our modular and data-driven approach for potent and selective covalent inhibitor discovery and provide a platform to apply it to other emerging targets.

Unmasking the Conformational Stability and Inhibitor Binding to SARS-CoV-2 Main Protease Active Site Mutants and Miniprecursor.

We recently demonstrated that inhibitor binding reorganizes the oxyanion loop of a monomeric catalytic domain of SARS CoV-2 main protease (MPro) from an unwound (E) to a wound (active, E*) conformation, independent of dimerization. Here we assess the effect of the flanking N-terminal residues, to imitate the MPro precursor prior to its autoprocessing, on conformational equilibria rendering stability and inhibitor binding. Thermal denaturation (Tm) of C145A mutant, unlike H41A, increases by 6.8 °C, relative to wild-type mature dimer. An inactivating H41A mutation to maintain a miniprecursor containing TSAVL[Q or E] of the flanking nsp4 sequence in an intact form [(-6)MProH41A and (-6*)MProH41A, respectively], and its corresponding mature MProH41A were systematically examined. While the H41A mutation exerts negligible effect on Tm and dimer dissociation constant (Kdimer) of MProH41A, relative to the wild type MPro, both miniprecursors show a 4-5 °C decrease in Tm and > 85-fold increase in Kdimer as compared to MProH41A. The Kd for the binding of the covalent inhibitor GC373 to (-6*)MProH41A increases ∼12-fold, relative to MProH41A, concomitant with its dimerization. While the inhibitor-free dimer exhibits a state in transit from E to E* with a conformational asymmetry of the protomers' oxyanion loops and helical domains, inhibitor binding restores the asymmetry to mature-like oxyanion loop conformations (E*) but not of the helical domains. Disorder of the terminal residues 1-2 and 302-306 observed in both structures suggest that N-terminal autoprocessing is tightly coupled to the E-E* equilibrium and stable dimer formation.

Autoprocessing and oxyanion loop reorganization upon GC373 and nirmatrelvir binding of monomeric SARS-CoV-2 main protease catalytic domain.

The monomeric catalytic domain (residues 1-199) of SARS-CoV-2 main protease (MPro1-199) fused to 25 amino acids of its flanking nsp4 region mediates its autoprocessing at the nsp4-MPro1-199 junction. We report the catalytic activity and the dissociation constants of MPro1-199 and its analogs with the covalent inhibitors GC373 and nirmatrelvir (NMV), and the estimated monomer-dimer equilibrium constants of these complexes. Mass spectrometry indicates the presence of the accumulated adduct of NMV bound to MProWT and MPro1-199 and not of GC373. A room temperature crystal structure reveals a native-like fold of the catalytic domain with an unwound oxyanion loop (E state). In contrast, the structure of a covalent complex of the catalytic domain-GC373 or NMV shows an oxyanion loop conformation (E* state) resembling the full-length mature dimer. These results suggest that the E-E* equilibrium modulates autoprocessing of the main protease when converting from a monomeric polyprotein precursor to the mature dimer.

EWALD: A macromolecular diffractometer for the second target station.

Revealing the positions of all the atoms in large macromolecules is powerful but only possible with neutron macromolecular crystallography (NMC). Neutrons provide a sensitive and gentle probe for the direct detection of protonation states at near-physiological temperatures and clean of artifacts caused by x rays or electrons. Currently, NMC use is restricted by the requirement for large crystal volumes even at state-of-the-art instruments such as the macromolecular neutron diffractometer at the Spallation Neutron Source. EWALD's design will break the crystal volume barrier and, thus, open the door for new types of experiments, the study of grand challenge systems, and the more routine use of NMC in biology. EWALD is a single crystal diffractometer capable of collecting data from macromolecular crystals on orders of magnitude smaller than what is currently feasible. The construction of EWALD at the Second Target Station will cause a revolution in NMC by enabling key discoveries in the biological, biomedical, and bioenergy sciences.

Cryo neutron crystallography demonstrates influence of RNA 2'-OH orientation on conformation, sugar pucker and water structure.

The ribose 2'-hydroxyl is the key chemical difference between RNA and DNA and primary source of their divergent structural and functional characteristics. Macromolecular X-ray diffraction experiments typically do not reveal the positions of hydrogen atoms. Thus, standard crystallography cannot determine 2'-OH orientation (H2'-C2'-O2'-HO2' torsion angle) and its potential roles in sculpting the RNA backbone and the expansive fold space. Here, we report the first neutron crystal structure of an RNA, the Escherichia coli rRNA Sarcin-Ricin Loop (SRL). 2'-OD orientations were established for all 27 residues and revealed O-D bonds pointing toward backbone (O3', 13 observations), nucleobase (11) or sugar (3). Most riboses in the SRL stem region show a 2'-OD backbone-orientation. GAGA-tetraloop riboses display a 2'-OD base-orientation. An atypical C2'-endo sugar pucker is strictly correlated with a 2'-OD sugar-orientation. Neutrons reveal the strong preference of the 2'-OH to donate in H-bonds and that 2'-OH orientation affects both backbone geometry and ribose pucker. We discuss 2'-OH and water molecule orientations in the SRL neutron structure and compare with results from a solution phase 10 µs MD simulation. We demonstrate that joint cryo-neutron/X-ray crystallography offers an all-in-one approach to determine the complete structural properties of RNA, i.e. geometry, conformation, protonation state and hydration structure.

The mechanisms of catalysis and ligand binding for the SARS-CoV-2 NSP3 macrodomain from neutron and x-ray diffraction at room temperature.

The nonstructural protein 3 (NSP3) macrodomain of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) (Mac1) removes adenosine diphosphate (ADP) ribosylation posttranslational modifications, playing a key role in the immune evasion capabilities of the virus responsible for the coronavirus disease 2019 pandemic. Here, we determined neutron and x-ray crystal structures of the SARS-CoV-2 NSP3 macrodomain using multiple crystal forms, temperatures, and pHs, across the apo and ADP-ribose-bound states. We characterize extensive solvation in the Mac1 active site and visualize how water networks reorganize upon binding of ADP-ribose and non-native ligands, inspiring strategies for displacing waters to increase the potency of Mac1 inhibitors. Determining the precise orientations of active site water molecules and the protonation states of key catalytic site residues by neutron crystallography suggests a catalytic mechanism for coronavirus macrodomains distinct from the substrate-assisted mechanism proposed for human MacroD2. These data provoke a reevaluation of macrodomain catalytic mechanisms and will guide the optimization of Mac1 inhibitors.

Hit Expansion of a Noncovalent SARS-CoV-2 Main Protease Inhibitor.

Inhibition of the SARS-CoV-2 main protease (Mpro) is a major focus of drug discovery efforts against COVID-19. Here we report a hit expansion of non-covalent inhibitors of Mpro. Starting from a recently discovered scaffold (The COVID Moonshot Consortium. Open Science Discovery of Oral Non-Covalent SARS-CoV-2 Main Protease Inhibitor Therapeutics. bioRxiv 2020.10.29.339317) represented by an isoquinoline series, we searched a database of over a billion compounds using a cheminformatics molecular fingerprinting approach. We identified and tested 48 compounds in enzyme inhibition assays, of which 21 exhibited inhibitory activity above 50% at 20 µM. Among these, four compounds with IC50 values around 1 µM were found. Interestingly, despite the large search space, the isoquinolone motif was conserved in each of these four strongest binders. Room-temperature X-ray structures of co-crystallized protein-inhibitor complexes were determined up to 1.9 Å resolution for two of these compounds as well as one of the stronger inhibitors in the original isoquinoline series, revealing essential interactions with the binding site and water molecules. Molecular dynamics simulations and quantum chemical calculations further elucidate the binding interactions as well as electrostatic effects on ligand binding. The results help explain the strength of this new non-covalent scaffold for Mpro inhibition and inform lead optimization efforts for this series, while demonstrating the effectiveness of a high-throughput computational approach to expanding a pharmacophore library.

Covalent narlaprevir- and boceprevir-derived hybrid inhibitors of SARS-CoV-2 main protease.

Emerging SARS-CoV-2 variants continue to threaten the effectiveness of COVID-19 vaccines, and small-molecule antivirals can provide an important therapeutic treatment option. The viral main protease (Mpro) is critical for virus replication and thus is considered an attractive drug target. We performed the design and characterization of three covalent hybrid inhibitors BBH-1, BBH-2 and NBH-2 created by splicing components of hepatitis C protease inhibitors boceprevir and narlaprevir, and known SARS-CoV-1 protease inhibitors. A joint X-ray/neutron structure of the Mpro/BBH-1 complex demonstrates that a Cys145 thiolate reaction with the inhibitor's keto-warhead creates a negatively charged oxyanion. Protonation states of the ionizable residues in the Mpro active site adapt to the inhibitor, which appears to be an intrinsic property of Mpro. Structural comparisons of the hybrid inhibitors with PF-07321332 reveal unconventional F···O interactions of PF-07321332 with Mpro which may explain its more favorable enthalpy of binding. BBH-1, BBH-2 and NBH-2 exhibit comparable antiviral properties in vitro relative to PF-07321332, making them good candidates for further design of improved antivirals.

Covalent narlaprevir- and boceprevir-derived hybrid inhibitors of SARS-CoV-2 main protease: room-temperature X-ray and neutron crystallography, binding thermodynamics, and antiviral activity.

The COVID-19 pandemic continues to disrupt everyday life, with constantly emerging SARS-CoV-2 variants threatening to render current vaccines ineffective. Small-molecule antivirals can provide an important therapeutic treatment option that is subject to challenges caused by the virus variants. The viral main protease (M pro ) is critical for the virus replication and thus is considered an attractive drug target for specific protease inhibitors. We performed the design and characterization of three reversible covalent hybrid inhibitors BBH-1, BBH-2 and NBH-2, whose structures were derived from those of hepatitis C protease inhibitors boceprevir and narlaprevir. A joint X-ray/neutron structure of the M pro /BBH-1 complex demonstrated that a Cys145 thiolate reaction with the inhibitor’s keto-warhead creates a negatively charged oxyanion, similar to that proposed for the M pro -catalyzed peptide bond hydrolysis. Protonation states of the ionizable residues in the M pro active site adapt to the inhibitor, which appears to be an intrinsic property of M pro . Structural comparisons of the hybrid inhibitors with PF-07321332 revealed unconventional interactions of PF-07321332 with M pro which may explain its more favorable enthalpy of binding and consequently higher potency. BBH-1, BBH-2 and NBH-2 demonstrated comparable antiviral properties in vitro relative to PF-07321332, making them good candidates for further design of improved antivirals.

The mechanisms of catalysis and ligand binding for the SARS-CoV-2 NSP3 macrodomain from neutron and X-ray diffraction at room temperature.

The NSP3 macrodomain of SARS CoV 2 (Mac1) removes ADP-ribosylation post-translational modifications, playing a key role in the immune evasion capabilities of the virus responsible for the COVID-19 pandemic. Here, we determined neutron and X-ray crystal structures of the SARS-CoV-2 NSP3 macrodomain using multiple crystal forms, temperatures, and pHs, across the apo and ADP-ribose-bound states. We characterize extensive solvation in the Mac1 active site, and visualize how water networks reorganize upon binding of ADP-ribose and non-native ligands, inspiring strategies for displacing waters to increase potency of Mac1 inhibitors. Determining the precise orientations of active site water molecules and the protonation states of key catalytic site residues by neutron crystallography suggests a catalytic mechanism for coronavirus macrodomains distinct from the substrate-assisted mechanism proposed for human MacroD2. These data provoke a re-evaluation of macrodomain catalytic mechanisms and will guide the optimization of Mac1 inhibitors.

Cryotrapping peroxide in the active site of human mitochondrial manganese superoxide dismutase crystals for neutron diffraction.

Structurally identifying the enzymatic intermediates of redox proteins has been elusive due to difficulty in resolving the H atoms involved in catalysis and the susceptibility of ligand complexes to photoreduction from X-rays. Cryotrapping ligands for neutron protein crystallography combines two powerful tools that offer the advantage of directly identifying hydrogen positions in redox-enzyme intermediates without radiolytic perturbation of metal-containing active sites. However, translating cryogenic techniques from X-ray to neutron crystallography is not straightforward due to the large crystal volumes and long data-collection times. Here, methods have been developed to visualize the evasive peroxo complex of manganese superoxide dismutase (MnSOD) so that all atoms, including H atoms, could be visualized. The subsequent cryocooling and ligand-trapping methods resulted in neutron data collection to 2.30 Šresolution. The P6122 crystal form of MnSOD is challenging because it has some of the largest unit-cell dimensions (a = b = 77.8, c = 236.8 Å) ever studied using high-resolution cryo-neutron crystallography. The resulting neutron diffraction data permitted the visualization of a dioxygen species bound to the MnSOD active-site metal that was indicative of successful cryotrapping.

High-Throughput Virtual Screening and Validation of a SARS-CoV-2 Main Protease Noncovalent Inhibitor.

Despite the recent availability of vaccines against the acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the search for inhibitory therapeutic agents has assumed importance especially in the context of emerging new viral variants. In this paper, we describe the discovery of a novel noncovalent small-molecule inhibitor, MCULE-5948770040, that binds to and inhibits the SARS-Cov-2 main protease (Mpro) by employing a scalable high-throughput virtual screening (HTVS) framework and a targeted compound library of over 6.5 million molecules that could be readily ordered and purchased. Our HTVS framework leverages the U.S. supercomputing infrastructure achieving nearly 91% resource utilization and nearly 126 million docking calculations per hour. Downstream biochemical assays validate this Mpro inhibitor with an inhibition constant (Ki) of 2.9 µM (95% CI 2.2, 4.0). Furthermore, using room-temperature X-ray crystallography, we show that MCULE-5948770040 binds to a cleft in the primary binding site of Mpro forming stable hydrogen bond and hydrophobic interactions. We then used multiple µs-time scale molecular dynamics (MD) simulations and machine learning (ML) techniques to elucidate how the bound ligand alters the conformational states accessed by Mpro, involving motions both proximal and distal to the binding site. Together, our results demonstrate how MCULE-5948770040 inhibits Mpro and offers a springboard for further therapeutic design.

Michaelis-like complex of SARS-CoV-2 main protease visualized by room-temperature X-ray crystallography.

SARS-CoV-2 emerged at the end of 2019 to cause an unprecedented pandemic of the deadly respiratory disease COVID-19 that continues to date. The viral main protease (Mpro) is essential for SARS-CoV-2 replication and is therefore an important drug target. Understanding the catalytic mechanism of Mpro, a cysteine protease with a catalytic site comprising the noncanonical Cys145-His41 dyad, can help in guiding drug design. Here, a 2.0 Šresolution room-temperature X-ray crystal structure is reported of a Michaelis-like complex of Mpro harboring a single inactivating mutation C145A bound to the octapeptide Ac-SAVLQSGF-CONH2 corresponding to the nsp4/nsp5 autocleavage site. The peptide substrate is unambiguously defined in subsites S5 to S3' by strong electron density. Superposition of the Michaelis-like complex with the neutron structure of substrate-free Mpro demonstrates that the catalytic site is inherently pre-organized for catalysis prior to substrate binding. Induced fit to the substrate is driven by P1 Gln binding in the predetermined subsite S1 and rearrangement of subsite S2 to accommodate P2 Leu. The Michaelis-like complex structure is ideal for in silico modeling of the SARS-CoV-2 Mpro catalytic mechanism.

Structural, Electronic, and Electrostatic Determinants for Inhibitor Binding to Subsites S1 and S2 in SARS-CoV-2 Main Protease.

Creating small-molecule antivirals specific for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) proteins is crucial to battle coronavirus disease 2019 (COVID-19). SARS-CoV-2 main protease (Mpro) is an established drug target for the design of protease inhibitors. We performed a structure-activity relationship (SAR) study of noncovalent compounds that bind in the enzyme's substrate-binding subsites S1 and S2, revealing structural, electronic, and electrostatic determinants of these sites. The study was guided by the X-ray/neutron structure of Mpro complexed with Mcule-5948770040 (compound 1), in which protonation states were directly visualized. Virtual reality-assisted structure analysis and small-molecule building were employed to generate analogues of 1. In vitro enzyme inhibition assays and room-temperature X-ray structures demonstrated the effect of chemical modifications on Mpro inhibition, showing that (1) maintaining correct geometry of an inhibitor's P1 group is essential to preserve the hydrogen bond with the protonated His163; (2) a positively charged linker is preferred; and (3) subsite S2 prefers nonbulky modestly electronegative groups.